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ABSTRACT 

The Pierre Auger Observatory reports that 20 of the 27 highest energy cosmic rays have arrival 
directions within 3.2° of a nearby galaxy in the Veron-Cetty & Veron Catalog of Quasars and Active 
Galactic Nuclei (12 th Ed.), with ~ 5 of the correlations expected by chance. In this paper we examine 
the correlated galaxies to gain insight into the possible UHECR sources. We find that 14 of the 21 
correlated VCV galaxies are AGNs and we determine their bolometric luminosities. The remaining 7 
are primarily star- forming galaxies. The bolometric luminosities of the correlated AGNs are all greater 
than 5 x lO^ergs" 1 , which may explain the absence of UHECRs from the Virgo region in spite of 
, the large number of VCV galaxies in Virgo, since most of the VCV galaxies in the Virgo region are 

^vq ■ low luminosity AGNs. Interestingly, the bolometric luminosities of most of the AGNs are significantly 

lower than required to satisfy the minimum condition for UHECR acceleration in a continuous jet. 
If a UHECR- AGN correlation is substantiated with further statistics, our results lend support to the 
i-jh ■ recently proposed "giant AGN flare" mechanism for UHECR acceleration. 

Subject headings: cosmic rays, AGN, X-rays 
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1. INTRODUCTION 



The origin of ultra- high en ergy cosmic rays (U HECRs) has been an important question in astrophysics, for many 
decades. Early observations (|Takeda et a l. 1998) suggested a violation of the Grizsen-Zatsepin-Kuzmin (GZK) pre- 
diction that the UHECR energy spectrum must drop off at energies above ~ 60 EeV. If confirmed, this would have 
been a major challenge for theory, but recent high-statistics observations by the^ HiRes and P ierre Auger coll abora- 
j^. | tions find a downturn in the spectrum consistent with the GZK prediction (|Abbasi et al J 12008 at lAbraham et alJ l2008) . 
t^J- , Nonetheless, the puzzle of which astrophysical sites are capable of accelerating UHECRs remains open, and the answer 
£C) • will be of fundamental importance for our understanding of gamma ray bursts (GRBs), active galactic nuclei (AGNs) 
and quasars, and other extreme systems. Efforts to find angular correlations between UHECR arrival directions and 
candidate astrophysical sources have been bedeviled until recently by a combination of inadequate statistics and the 
fact that UHECRs, being charged particles, are deflected by magnetic fields en route from their sources. 

The Pierre Auger Observatory's discovery of a significant correlation between the highest energy cosmic rays and 
nearby galaxies in the[V eroii-Ccttv & Veron (2006) (VCV) Catalog of Quasars and Active Galactic Nuclei (12 th Ed.) 
(The Pierre Auger Collabora tion 2007a. b) (Auger07a,b below) is an important step toward identifying the sources of 
UHECRs. Of the 27 cosmic rays above 57 EeV recorded prior to Aug. 31, 2007, twenty are within 3.2° of a VCV 
galaxy with z < 0.018 (about 75 Mpc). Restricting to \b\ > 10°, where the VCV catalog is more complete, there are 22 
UHECRs of which 19 are correlated. About 5 of these correlations would be expected by chance if the arrival directions 
of UHECRs were isotropic. The distribution of magnetic deflections is not known a priori and it is not yet known with 
certainty whether the UHECRs are protons or nuclei. Consequently the expected correlation between UHECR energy 
and maximum source distance is not certain. The strategy used by Auger was to search for correlations by scanning 
over UHECR energy threshold, maximum angular separation, and maximum source redshift, to find the values of these 
parameters that maximize the significance of the correlations. The VCV catalog and the parameters given above were 
identified using data through May 31, 2006 and a "prescription" formulated for an independent test. The correlation 
was confirmed with independent subsequent data, taken from June 1, 2006 to Aug. 31, 2007, with a probability of less 
than 1% of occurring by chance. 

This Auger result is of fundamental importance to particle astrophysics, because the correlation with nearby extra- 
galactic structure clearly demonstrates that UHECRs are of extragalactic origin and that the highest energy cosmic 
rays have a horizon, consistent with the GZK effect. (The downturn in the spectrum might merely be due to a max- 
imum energy of the accelerators.) However as stressed by the Auger collaboration, the observed correlation may not 
mean that the correlated UHECRs are produced by galaxies with which they are correlated: the VCV galaxies may 
just be tracers of the true sources. 

Our purpose here is to examine the VCV galaxies that correlate with the 20 Auger UHECRs, as a step toward 
elucidating properties of the sources of UHECRs. We emphasize that the correlation observed by Auger is only 
statistical: ~l/4 of the correlations are expected simply by chance. We cannot be confident that any given one of 
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the correlated VCV galaxies is a source. However the degree of correlation observed by Auger is so high that galaxy 
clustering alone cannot account for it (GRF, A. Bcrlind and IZ, in preparation); of order half or more of the correlated 
galaxies are most likely the sources of their associated UHECR, enabling us to obtain statistically useful information 
on source properties even though we do not have a pure sample of sources. The underlying principle is that the chance 
correlations must be representative of the ensemble of VCV galaxies, and therefore their presence does not distort the 
conclusions when properly used. 

The VCV catalog is a list of the active galactic nucleus, quasar and BL Lac candidates reported in the literature, 
based on heterogeneous selection criteria. While it is the largest available collection of known AGNs, especially for the 
southern hemisphere where most of Auger's exposure lies, it has several deficiencies compared to an ideal catalog. The 
VCV catalog is known to be incomplete and non-uniform. Furthermore, it is not pure. Here we assemble the correct 
classifications of the VCV galaxies correlated with the Auger UHECRs. In addition, we determine the bolometric 
luminosity of each AGN, which is an important diagnostic of the maximum UHECR energy an AGN can produce 
in conventional jet acceleration. Finally, we look for evidence of a bolometric luminosity threshold for AGNs which 
may be responsible for UHECR production. We examine the AGNs in the Virgo region, from which no UHECRs 
were detected, and find that most are much lower in their bolometric luminosity than the AGNs which are correlated 
with UHECRs. Thus the lack of UHECRs from the direction of Virgo may simply reflect the existence of a minimum 
bolometric luminosity of AGNs responsible for accelerating UHECRs. We also estimate the fraction of low-luminosity 
AGNs (LLAGNs, taken here to have Lboi < 5 x 10 42 ergs _1 ) in the VCV catalog in the HiRes exposure region, finding 
that about half the AGNs are below the luminosity threshold of AGNs correlated with Auger UHECRs. The fraction of 
known L LAGNs in VCV is higher in the HiRes field of view than in Auger's, due to the sensitive northern hemisphere 
survey by I Ho et al.l (|199E ); this may contribu te to HiRes' not observing a significant correlation between VCV galaxies 



and their UHECR data (|Abbasi et al.ll2008bt) 

2. UHECR ACCELERATION 

The motivation and agenda for our study derive from theoretical considerations of UHECR acceleration. No as- 
trophysical system has been conclusively demonstrated theoretically to be capable of ac celerating cosmic ray s to the 
observed energies of > 10 20 eV. GRBs have been argued to be responsible for UHECRs (jWaxmaril 19951 l2004f ) on the 
basis of their energy injection rate and theoretical plausibility: they are known to produce high energy photons and the 
GRB internal-shock model can be viable for UHECR acceleration as well. Other possible accelerators include internal 
shocks in the jets of Active Galactic Nuclei - analogous to those in GRBs but with much lower bulk Lorentz factors 
(for an early su ggestion see 



(tor an early su ggestion see i Biermann fc Strittmatterl (|1987| )h external shocks in the lobes of powerful radio galaxies, 
and magnetars (jAronsI 12003 ), to name s ome of the more popular. A GRB could satisfy the requirements of UHECR 



acceleration for T ~ 10 3 ()Waxmanlfl995T l. while for AGN jet acceleration T ~ few would be envisaged. Assuming the 
correlation between UHECRs and VCV galaxies observed by Auger is real, AGNs become the favored sources and we 
focus here on testing AGN-based acceleration models. 

A very general requirement to accelerate a cosmic ray proton of energy E = £20 10 20 eV, in a re lativistic jet of bulk 
Lorentz factor T, is that the Poynting luminosity of the jet, Lp, satisfy ip > 10 45 T 2 E 2 ergs -1 (jFarrar fc Gruzinovl 
2008). This follows because the Larmour radius of the UHECR 1 , 0.1 Mpc i?2o/^G: must be less than the characteristic 
size of the acceleration region, R, leading to RB>3 x 10 17 T" 1 E 20 G cm, which translates to a lower bound on the 
Poynting luminosity of the jet, Lp ~ \cT A B 2 R 2 . For a conventional AGN jet the Poynting luminosity is supplied by 
accretion, so both the jet an d accretion disk of t he ac tive nucleus must satisfy the quoted luminosity constraint, as 
discussed in greater detail bv lFarrar fc Gruzinovl ((2008) . Therefore, the bolometric luminosity of the accretion disk of 
an AGN capable of accelerating a proton to E20 10 20 eV must satisfy 

L bol > 10 45 r 2 L 2 ergs- 1 . (1) 

If the time-scale of variation of the accele rator is large compared to the cosmic ray travel time, estimated to be < 10 5 
yr for the correlated Auger cosmic rays (jFarrar fc Gruzinovl [2008) . then the required luminosity should be evident 
when observing the sources today. In conventional AGN jet acceleration or acceleration in radio lobes, the time scale 
for variation of the source is of order the lifetime of an AGN, i.e., > 10 7 yr. Thus an AGN capable of accelerating 
UHECRs should either have powerful radio lobes or have an accretion disk luminosity consistent with the bound on 
Z/boi in equation JTJ). 

An alternative t o the conventional continuous jet model of AGN acceleration has recently been proposed by 
iFarrar fc Gruzinovl (|2008f) . In this "Giant AGN Flare" mechanism, an instability of the accretion disk - perhaps 
initiated by the tidal disruption of a passing star - produces an intense flare l asting of order a day to a month. During 
the flare the luminosity condition for UHECR acceleration is easily satisfied (jFarrar fc Gruzino v 2008) , but afterward 
the emission subsides quickly because the cooling time is short and the material which fuels the flare is largely con- 
sumed, so the system as observed today need not satisfy the luminosity condition |T]) . In section [5] we determine, for 
each AGN-UHECR pair, the figure-of-merit Aboi = LboilO -45 -^ 2 - Aboi should be > 1 for conventional continuous 
jet acceleration of protons, but has no such constraint for the Giant Flare scenario. We find that at least half of the 
correlated AGNs do not satisfy the bolometric luminosity requirement, favoring the Giant AGN Flare scenario. 

A question raised (but not yet answered) by the Giant AGN flare model, is whether the AGNs responsible for 
UHECR production, as observed today, show a threshold Lboi- In the AGN flare model an accretion disk is required 
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prior to the flare, but the observed value of Lboi reflects the properties of the remnant system, which depends on how 
much material is left after the period of rapid accretion, and how cool the accretion disk has become in the intervening 
time. In section [S] we find empirical evidence of a threshold Lboi > 5 x 10 42 ergs . 
We now turn to our task; our first objective is to identify the AGNs among the correlated VCV galaxies. 

3. IDENTIFICATION OF ACTIVE GALACTIC NUCLEI 

There are a number of diagnostics used to identify an accreting supermassive Black Hole from optical spectroscopy. 
Broad permitted emission lines (e.g. Ha) with line widths > 1000 km/s strongly suggest that we are seeing the 
Keplerian velocities of gas orbiting the BH. This dense orbiting gas, known as the broad-line region, has a typical size 
of light days, and thus is never spatially resolved in astronomical observations. Objects with observable broad emission 
lines, for which we have an unobstructed view of the nucleus, are known as broad-line AGNs (or, at low luminosities, 
Seyfert 1 galaxies). In some cases our view of the central engine and broad-line region are obscured. AGNs may still 
be identified by the presence of narrow emission lines, originating from more rarefied gas at scales of 100-1000s of pc. 
In particular, the ratios of strong emission lines (typically H/3, [O III] A5007, Ha, and [N II] A6584) are very sensitive 
to the shape of the ionizing continuum, and thus allow us to distinguish between nebulae excited by starlight, shocks, 
or the hard continuum of a radiating BH. Obscured objects are known as narrow-line AGNs (Seyfert 2 galaxies at low 
luminosities). 

Two-dimensional line ratio diagnostics have been developed to effectively discrim inate between variou s emission 
mecha nisms, e.g., star formation, s hock, or photoionization by an accreting B H. See iBaldwin et all ([1981) (BPT be- 
low), (jVeilleux k Osterbrocklll987l : iKewlev et all I2TM I Kauffmann et al.ll2003h . Baldwin, Phillips, k Terlevic (BPT) 
diagnostic diagrams are based on the relative strengths of prominent emission lines (e.g. [O III] /H/3 versus [N II] /Ha) 
that are close together in wavelength space, to minimize the impact of reddening. A typical BPT diagra m is shown in 
Figur e [TJ along with various boundary lines between different types of galaxies with narrow-line emission. IKewlev et al.l 
(|2001[ ) developed a conservative boundary between star formation and AGN activity based on theoretical modeling, 
shown as a solid curve in Figure [TJ galaxi es lying above and to the right of the Kewley line are unambiguously AGNs. 
A looser, empirical boundary was used by Kauffm ann et al.l (|2003l ) to identify SDSS star-forming galaxies by their line 
ratios; it is shown as a dashed curve in Figure [TJ There is a further distinction for the objects which fall outside the 
star formation boundary: those with log([0 III] /H/3) > 0.48, are known as high-ionization Seyf ert galaxies, and those 
with log([0 III]/H/3) < 0.48 are known as low-ionization nuclear emission regions (LINERs) (|Heckmanlll980l) . The 
Seyfert-LINER boundary is indicated by the horiz ontal line in Figure [TJ While many LINERs have been demonstrated 
to be powered by accretion activity (see review in I Hoi (|2008h ). there are other processes such as shocks (Veilleux et al. 
[19951) which can lead to LINER-like line ratios. Caution, and auxiliary data, are often required to determine the 
nature of distant LINERs. 

It is, of course, important to remember that all AGN selection techniques are biased in different ways, and are 
dependent on the depths and apertures of in dividual surveys. For nearby AGNs, optical emiss ion line surveys are 
proba bly the most complete, e.g., as shown by (|Heckman et al.|[2005l : Ho 2008). Note, however, that lReviglio k Helfandl 
(2006) find that about half the AGNs detected in their sample using radio or X-ray selection are not identified by the 
BPT criteria discussed above. Since UHECR source candidates are closer than about 100 Mpc, and star-forming regions 
with enough emission to be confused with an AGN are larger than « 10 pc, X-ray observation with Chandra resolution 
can decide whether a UHECR candidate source is an AGN or not. In the radio, detection of jets unambiguously 
identifies an AGN. If no jets are seen, the presence of an AGN can be inferred if there is a compact source at the 
center of the galaxy whose radio emission excee ds that which would be expected for a nuclear starburst region based 
on the far infrared-radio relation (|Condonlll992f ). 

4. UHECR-CORRELATED VCV GALAXIES 

Auger has published 27 UHECRs with energies above 57 EeV, 20 of which are correlated with a VCV galaxy closer 
than z — 0.018 (Auger07a,b). That is our basic UHECR sample for this study. The 20 UHECRs are correlated with 
21 VCV galaxies. 2 Table [TJ shows the 20 correlated UHECRs and the 21 associated VCV galaxies. For each VCV 
galaxy, we check the VCV classification against literature and by examining available spectra. Five are unambiguously 
broad-line AGNs and one is a BL Lac - widely considered to be a radio-loud AGN viewed directly down the jet. 
The remaining 15 galaxies are potential narrow- line AGNs and have to be evaluated individually. Figure [TJ is a 
BPT diagnostic diagram in which the line ratios, log([OIII]/H/3) vs. log([NII]/Ha), for 11 of the 15 galaxies are 
plotted. The four galaxies which cannot be plotted in Fig. [T] consist of two AGNs and two without nuclear activity. 
The optical emission from NGC 4945 is highly obscured but its hard X-r ay spectrum clearly i dentifies AGN activity 
dMadeiski et al.ll2000h. E SQ 139-G12 sho ws a hint of a broad ened Ha line ([Marquez et al.ll2004h and is an AGN. NGC 
5244 (M oran et al.lll996l ) and NGC 7135 (jJoguet et al.ll2001T ) have almost no [OIII] emission; they cannot be identified 
as AGNs from their optical spectra. 

Eleven galaxies appear in Figure [TJ of which six, shown as diamonds on the BPT diagram, fall comfortably in the 
Seyfert region of the BPT diagram; these a re clear AGNs. Two of the eleven (NGC 7591 and Q 2207+0122), shown 
as a triangle and an asterisk, fall below the [Kauffmann et al. (2003) empirical star formation line and are not AGNs. 
The remaining three (NGC 2989, NGC 1204, I C 5169), shown as tr iangles, fall between the Kewley and Kauffmann 
lines, in the LINER region of the BPT diagram. IKewlev et all ([2001 uses additional information from other emission 
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TABLE 1 

Auger UHECRs and Correlated VCV Galaxies 
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Note. — This table lists each Auger UHECR (with E > 57 EeV) whose arrival directions are within 3.2° of a nearby (z < 0.018) VCV 
galaxy, along with the correlated VCV object. The year and the Julian day when the UHECRs were recorded are given, as well their 
energies and positions (equatorial and galactic) in degrees, from Auger07b. This is followed by the name of the VCV galaxy, its position 
in degrees, and its separation from the UHECR in degrees, its rcdshift, and the VCV classification." The last column shows the correct 
classification (taken from the literature where available and confirmed in every case). The UHECRs correlated with galaxies that are 
shown here to have active nuclei are marked with an asterisk. 

Ref ere nces. — (1) ISchneider et al I JT994T) , (2) jKewlev et al.l (f200lT). ( 3) ISturm et al.l <200(3). IVeilleux efahl tl99Sfi . (4) IMoran et all 
1T99H) , (5) IConcalves et al.l 119991) , (6) ISturm et al.l (120061 ), IVeilleux et all 1119951 ), (7) IJoeuet et al.l [1200111 

"Star forming regions are also known as HII or H2 regions because young stars ionize the hydrogen around them (neutral hydrogen is 
known as HI). 

lines and determines IC 5169 to be purely star forming; NGC 2989 ([Phillips et al.lll983f) and NGC 1204 (jCorbett et all 
2003) are predominantly star-forming galaxies but may have limited AGN activity, contributing < 25% of the total 
luminosity. Detailed information on each of the correlated VCV galaxies is given in Appendix \X\ and summarized in 
Tables [2] and [3 

In all, 14 of the 21 VCV galaxies, correlated with 14 UHECRs, are AGNs, while seven show limited or no sign 
of AGN activity. Radio or X-ray observations of these 7 VCV galaxies which do not appear to be AGN are needed 
to determine whether or not all UHECRs may be attributable to AGNs. 3 Figure [2] shows Auger's 27 highest-energy 
cosmic rays with the correlated UHECRs represented by 3.2° circles and the uncorrelated ones represented by 8-point 
stars. The 14 confirmed AGNs are shown as triangles and the 7 other correlating VCV galaxies are shown as squares. 
Auger's exposure cutoff at S — 25° is indicated (cyan line) as is the Galactic plane band = ±10°), where VCV is 
especially incomplete. The rest of the VCV galaxies with z < 0.018 are shown as points, color-coded in bins of redshift. 

5. BOLOMETRIC LUMINOSITIES 

Under the widely-used assumption of a universal spectral energy distribution for active galaxies, the bolometric 
luminosities may be estimated from measurements of a single line or in a single band. In most cases, we have 
flux measurements in only one waveband. In this paper, we use the following bolometr ic corrections for [011 1] and 
broad Ha emission: iboi = 3500 £[oiii]j with a 0.38 dex variance, as determined by iHeckman et ahl |2004), and 
L bo i = {L(Ha broad )/10 42 )°- S6 x 2.34 x lO^ergs" 1 dGreene fc Hoi [2007m . 

The conversion fa c tor fr om the 2-10 keV X-ray luminosity to the bolometric luminosity, as determined by 
iVasudevan fc Fabianl ([2007), is dependent upon the level of AGN activity; it is 15-25 for AGNs with Eddington 
ratios [Lboil ^Edd) below ~0.1, and 40-70 for those with Eddington ratios above ~0.1. The Eddington luminosity, 

3 A proposal has been submitted to Chandra by GRF, W. Terrano and IZ to make the observations required. 
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Fig. 1.— A BPT diagram showing the line ratios, [OIII]/H/3 vs. [NII]/Ha, for 11 of the 15 VCV galaxies correlated with UHECRs which 
are not a broad-line AGN or BL Lac. NGC 5244 and NGC 7135 are not plotted since they show almost no [OIII] emission, nor are NGC 
4945 and ESO 139-G12 since their identification as AGNs is made on other grounds. The dwarf galaxy, Q 2207+0122, is indicated by the 
asterisk. NGC 2989, IC 516 9, NGC 1204, and N GC 7591 are shown as triangles. These fall to the left and below the theoretical classification 
line for AGN developed by Kcwlcy ct al. (2001), indicated by the solid li ne. The Seyfert 2 galaxi es are shown as diamonds. They lie well 
within the Seyfert region of the diagram. The demarcation line used by Kauffmann ct al. ( 2003) is shown (dashed) for comparison. The 
horizontal line shows the boundary between Seyferts (above the line) and LINERs (below the line). 

Fig. 2. — Sky plot in equatorial coordinates showing Auger's 27 highest energy cosmic rays and VCV galaxies. The UHECRs correlating 
with an AGN are indicated with open 3.2° circles and those that do not are indicated with an 8-point star. The 14 correlated VCV galaxies 
which are known to be AGNs are shown as triangles. The 7 other correlating VCV galaxies are shown as squares. The rest of the VCV 
galaxies with z < 0.018 are shown as points color-coded by redshift. The Galactic plane band (|6| = 10°) where the VCV catalog is most 
incomplete is indicated by the solid black lines. Auger's exposure cutoff at <5 < 25° is indicated by the cyan line. 

LEdd, is related to the SMBH mass, M, by the relation L Edd = 3.3 x 1O 4 (M/M )L . We were able to obtain SMBH 
mass estimates for four AGNs in ou r sample; two can be found in the literature and we derive two from the broad Ha 
luminosity and line width following I Greene &: Hoi (|2005l ). (See details in Appendix A.) All four have low Eddington 
ratios: 0.006, 0.015, 0.029, and 0.1. Therefore, in the absence of SMBH mass measurements for the rest, we use 
Lboi = 20i 2 -iokoV for all. 

We determine the bolometric luminosities for the correlated AGNs from [OIII], broad Ha, or 2-10 keV emission. For 
the broad line AGNs, the [OIII] and Ha luminosities are taken from published data when possible, or fit from available 
optical spectra. (See Appendix A for details.) We use multicomponent Gaussian models to fit for the broad Ha 
luminosities. When fitting for [OIII] luminosities, we assume only a single Gaussian. The 2-10 keV X-ray luminosities 
are collected from literature. There is one AGN (ESO 139-G12), for which neither a calibrated spectrum nor an X-ray 
detection is available, and we use the upper limit for the X-ray luminosity. The AGNs and their bolometric luminosities 
are given in Table [2] 

For each AGN-UHECR pair, the figure-of-merit A bo i = LboilO -45 ^ 2 is determined. Out of the 13 UHECRs which 
are correlated with AGNs for which we were able to determine the bolometric luminosities, we see that only one 
UHECR correlates with an AGN having a value of Aboi ^ 1. Another 5 UHECRs are correlated with AGNs whose 
luminosity is within a factor of 2 of the minimum. Given the rough nature of the Aboi bound, these AGN may 
have sufficient power to accelerate their associated UHECRs. However seven correlated UHECRs are associated with 
AGNs that fall far short of the minimum power. This poses a serious problem for conventional models of UHECR 
acceleration in con t inuou s AGNs, although it is compatible with the recently proposed giant AGN flare scenario of 
iFarrar fc Gruzinovi (|2008f >. 

6. POSSIBLE THRESHOLD IN BOLOMETRIC LUMINOSITY 

We have seen in section [5] that only one of the AGNs correlated with UHECRs satisfies the naive limit for UHECR 
acceleration. We seek in this section to determine whether correlated AGNs show any threshold bolometric luminosity 
- such a threshold would elucidate the nature of the UHECR acceleration mechanism, and as we shall argue, can 
account for puzzling features of present observations. 

One puzzling aspect of the observed correlation of Auger UHECRs with nearby VCV galaxies is the lack of UHECRs 
from the Virgo Cluster (Auger07b). Because of its proximity, and the very large number of VCV galaxies in Virgo, 
this seems quite surprising. Therefore, we have made the same analysis as described above, for VCV galaxies in Virgo. 
We examined the V CV galaxies within a circle of radius 15° centered on the Virgo Cluster at aj2ooo = 12:26:32.1, 
<>J20oo = +12:43:24 (iBinggeli et al.lfl987f l and z < 0.009, which are attributed to the Virgo Cluster by (or are not 
listed in) iHo et all ([1997a|). There are a total of 30 VCV galaxies which satisfy these criteria. VCV classifies five 
of them as HII galaxies. One other galaxy was also found to be an HII galaxy. We can determine the bolometric 
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Bolometric Luminosities of Correlated AGN 
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NGC 4945 
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Note. — This table lists key properties of the VCV AGNs which are within 3.2° of an Auger UHECR,. The year, Julian day, and energy 
(in EeV) of the correlated UHECR are given for each AGN. The separation between the AGN and UHECR are given in degrees. The 
redshift and type of the AGN are given. (NLS1 indicates the presence of both narrow and broad line emission.) The derived infrared 
and bolometric luminosities are listed. \\,ol ls defined as Lj, ; X 10 -45 X E20 — 2 , where E20 = -E/100; it is the ratio which compares the 
bolometric luminosity to the theoretical luminosity necessary for producing the associated UHECR. 



luminosities of 23 out of th e remaining 24. In 18 cases we use the measured 2-10 keV X-ray lum i nosities (iHo et all 
20011 iTerashima et al.ll2002t ISatvapal etall 120051 : iGonzalez-Martm et al.|[2006t iPanessa et aLll2006t iHorst et alj 120081 : 
Kanda lvan 20 051), and for the five f or which X-ray luminosities are not available we use the measured [OIII] flux 
(|Shields et all 120071 lHolFaIlll997ah . There is some eviden c e tha t the [OIII] to bolometric luminosity conversion 
depends on luminosity or Eddington ratio (e.g., iNetzer et all (|2006[)'). S ince optical emission line fluxes, particularly 
Ha, may be contaminated by non-nuclear sources (e.g.. IHo et alT Eool) we have chosen to use the hard X-ray flux 
as a bolometric indicator whenever possible. 

The bolometric luminosities of the Virgo AGNs, listed in Appendix [B] are systematically lower than those of the 
AGNs correlated with UHECRs (which span the range 5 x 10 42 to 1 x 10 46 erg s" 1 ). Of the Virgo AGNs, only two 
have bolometric luminosities in the range of the correlated AGNs (NGC 4338 and NGC 3976, with luminosities of 1.1 
x 10 43 erg s" 1 and 4.6 x 10 42 respectively) and three others (NGC 4486, NGC 4579, and NGC 4698) have bolometric 
luminosities within a factor of five lower. The remaining 18 have bolometric luminosities ranging from 10 39 erg s _1 to 
10 41 erg s _1 , far below the bolometric luminosities of AGNs correlated with UHECRs. A histogram of the bolometric 
luminosities of the Virgo AGNs, compared with a histogram of the bolometric luminosities of the AGNs correlated 
with UHECRs, is shown in Figure [31 The abundance of known low luminosity AGNs in Virgo i s due to the sensitive 
Palomar spectroscopic survey of near by galaxies in the northern hemisphere bv lHo et alj (|1995l ). which found AGNs 
in ~60% of bulge-dominated galaxies (Ho "et~aLlll997bD . 

It may be that the accretion activity of these very low luminosity AGNs, is too low for them to be sources of UHECRs. 
A lower bound on the quiescent accretion rate might be expected if UHECRs are accelerated in intense, month-long, 
giant AGN flares, tr iggered when the tidal stre am of a star disrupted by a SMBH interacts with an existing accretion 
disk, as proposed bv lFarrar fc Gruzinovl (|2008l ). Note that if there is a threshold accretion luminosity, it may be larger 
than iboi.min of the correlated AGNs, since ~ 3 of the correlated AGN s are presumably chance associations, and the 
AGN luminosity function increases at low luminosities (|Ho et al.lll995l) . 

Another puz zle is the lack of sig nificant correlation between UHECRs and VCV galaxies reported by the HiRes 
Collaboration (lAbbasi et alj|2008b ). If UHECRs are preferentially accelerated by AGNs above a threshold bolometric 
luminosity, then one would expect a reduced degree of correlation significance observed by northern hemisphere exper- 
iments such as AGASA and HiRes, compared to what is seen by Auger in t he southe r n hem isphere. This is because 
a substantial fraction of the VCV northern hemisphere AGNs are LLAGNs IHo et all (fl995l) and 14% of Au ger's ex- 
posure and 84% of HiRes' exposure are from the northern hemisphere. We can estimate the degree of diluti on in the 
VCV catalog for z < 0.018 in the HiRes exposure compared to that in the Auger exposure, as follows. The lHo et al.1 
(1995) survey included 417 emission line galaxies with z < 0.018. They report that about 50% are AGN, with median 
L(Ha) = 2 x 10 39 erg s _1 , which roughly corresponds to our empirical UHECR threshold Lboi.mm = 5 x 10 erg 
s _1 . Thus at least « 105 of the VCV galaxies in the northern hemisphere with z < 0.018 are LLAGNs. The exposure 
weighted average fraction of LLAGNs and HII galaxies compared to total VCV galaxies is « 0.443 for HiRes, but only 
* 0.193 for Auger. This is a concrete illustration of the systematic differences in completeness and purity, between 
the portions of the VCV catalog used by the two UHECR experiments. 
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Fig. 3. — Histogram of the bolometric luminosities of Virgo AGNs (thin outline) compared with the bolometric luminosities of the AGNs 
correlated with UHECRs (thick outline). The correlated AGNs are much more luminous than the Virgo AGNs on average. 



TABLE 3 

Correlated VCV Galaxies with Limited AGN Activity 



AGN 


Auger Yr, Day 


E 
(EeV) 


r 

(deg) 


z 


log([NII]/H«) 


log([OIII]/H/3) 


log(L IR (8-1000 /Ltm)) 
(Lq) 


Type 


NGC 2989 


2007,84 


64 


3.00 


0.013 


-0.28 


0.23 


10.48 


Starburst (AGN < 20%) 


NGC 7591 


2006,185 


83 


3.04 


0.017 


-0.32 


-0.35 


11.07 


LIRG 


NGC 1204 


2007,145 


78 


1.57 


0.014 


-0.11 


-0.47 


10.83 


Starburst (AGN < 25%) 


NGC 5244 


2006,299 


69 


2.70 


0.008 






10.41 


Starburst 


IC 5169 


2005,295 


57 


2.13 


0.010 


-0.17 


0.20 


10.36 


Starburst 


NGC 7135 


2007,193 


90 


2.12 


0.007 






9.42 


Quiescent 


Q 2207+0122 


2005,63 
2007,51 


71 
58 


3.16 
1.58 


0.013 


-1.58 


0.82 




Dwarf Galaxy 



Note. — This table lists the VCV galaxies within 3.2° of an Auger UHECR that do not show evidence in their optical spectra of significant 
nuclear activity. The year, Julian day, and energy of the correlated UHECR as well as the redshift of the galaxy and its separation from the 
correlated UHECR are given. These are followed by line ratios ([Nil] /Ho and [OIII]/H/3), fit or gathered from literature. The line ratios 
for NGC 5244 and NGC 7135 are not available but b oth show almost no [0111] emission. The IR luminosities (8-1000 /im) are calculated 
from the IRAS fluxes according to the prescription of Sanders & Mirabel (1996). The final column classifies the galaxies according to their 
IR luminosities, which are indicators of their star formation rates. 

References. — NGC 2989 : IGoncalves etaLI (IT999T ). IC 5169: IKewlev et all (I2001T) , NGC 1204: I Veilleux et all (IT995T) , NGC 7591: 
Moustakas & Kcnnicutt (2006). The line ratios (courtesy of J. Moustakas) for Q 2207+0122 are fit from the spectrum (courtesy of J. 
Gunn and D. Schneider). 

7. INFRARED LUMINOSITIES 

An interesting question is whether the correlated VCV galaxies which are not optical AGNs have any common 
properties. One relevant property could be the total infrared luminosity, which is sensitive to both star formation and 
nuclear activity. In fact, there is circumstantial evidence that accretion is accompanied by significant star formation 
(e.g., Sanders & Mirabel 1996; Kauffmann et al. 2003), although the temporal coincidence of the two is a matter of 
debate (Ho 2005). It is certainly theoretically plausible that the large gas accretion episodes leading to AGN activity 
are accompanied by significant star formation events, and correspondingly elevated dust levels. Sufficiently dense and 
dusty star formation, such as is seen in ultra-luminous infrared galaxies, may in principle completely obscure AGN 
activity. 

We determi ne the IR luminosities of t he seven VCV galaxies that are not AGNs from IRAS fluxes following the 
prescription of ISanders fc Mirabel! (|1996[ k 

F m = 1.8 x 10- 14 [13.48/ 12 + 5.16/ 2 5 + 2.58/ 60 + /ioo] W/m 2 , L m = AnD 2 L F m , (2) 

where Dl is the luminosity distance of the galaxy, and / 12 to / 10 o are the IRAS fluxes. One of these non-optical- AGN 
galaxies is a luminous infrared galaxy (LIRG) (high rate of star formation: 11 < \og\Lmj Lq] < 12), four are starburst 
galaxies (10 < \og[Lm/ Lq\ < 11), one is a quiescent galaxy (\og{Lin/ L & ] < 10), and one has very low luminosity 
- it is in fact a dwarf galaxy. The infrared luminosities and classifications of these galaxies ar e listed in Table [3l 
For co mpleteness, we also calculate the IR luminosities of the AGN and include them in Table [5] ISanders'fe Mirabel 
(|1996h report that about 15% of LIRGs are Seyferts and another 35% are LINERS, whereas for starburst galaxies the 
LINER fraction is about the same but <5% are Seyferts. 



8. NOTE ADDED - IMoskalenko et al.l j2008) 



While this paper was being reviewed by the Pierre Auger Collaboration Pub l icatio n Committee (submitted on 
April 27, approved on May 28, 2008) a paper has appeared bv iMoskalenko et al.l (|2008l ) (MSPC below) which takes 
a complementary approach to ours, examining the radio morphology and luminosity of galaxies that might be the 
sources of the Auger UHECRs. MSPC extends the set of candidate sources to include 27 galaxies: 19 of the 21 VCV 
galaxies considered here, plus 8 more galaxies with z < 0.018 that are also candidate AGNs according to the NASA 
Extragalactic Database (NED) but which do not appear in the VCV catalog, plus an additional 27 galaxies from VCV 
and NED in the redshift range 0.018-0.037. The 19 VCV z < 0.018 galaxies that MSPC consider are those in our Table 
1, less SDSS J03349-0548 and ESO 383-G18. The former is removed because they exclude SDSS galaxies (except that 
SDSS J053302-0532 is in their sample as NGC 1346) but we do not understand why they do not include ESO 383-G18. 

Of the 54 galaxies considered, MSPC reports that 4 of them display the powerful jets or radio lobes expected for 
UHECR sources: 2 in each of their redshift bins. One of the two in the lower redshift bin, Cen A (NGC 5128) with 
Aboi = 0.03, is in our sample and the other, PKS 1343-60, is not. We examine the optical properties of PKS 1343-60 
found in literature. Although it has a bright optical nucleus, the emission is highly absorbed and it could not be 
determined as an AGN using opical criteria (See Appendix C for details). This most likely explains its absence from 
the VCV catalog which is largely composed of AGNs identified by optical criteria. We note that while Cen A has the 
largest radio flux of their candidate sources, 54 Jy, its extreme proximity makes it a weaker radio source than PKS 
1343-60 with 23 Jy. 

For completeness, we have extended the study reported in section [4] to include the 8 new AGN candidates with z < 
0.018 within 3.2° of an Auger UHECR introduced by MSPC. Only two have published optical line fluxes, WKK 4374 
(IGR J14515-5542) and NGC 7626. The former, WKK 4374, is a narrow line Seyfert (S2) with L bo i = 1.9 x lO^ergs" 1 
and Aboi = 0.46; intriguingly, it correlates with a previously uncorrelated UHECR (Auger year, day = 2007, 186; E=64 
EeV). The latter is a radio-loud LINER with L bo] = 1.2 x lO^ergs -1 , A bo i = 0.002. Two of the UHECRs whose 
VCV-correlated galaxies were found here to not be identifiable AGNs, (2006, 185; 83 EeV) and (2007, 84; 64 EeV), 
correlate with the new "MSPC" source candidates NGC 7626 and NGC 2907 respectively, although a spectrum is 
needed to determine whether the latter meets optical AGN criteria. 

The main result of MSPC is that apart from the four galaxies with powerful jets or radio lobes, their other source 
candidates have no special features that distinguish them from generic Seyferts and LINERs. On that basis, they 
conclude that the UHECR- AGN correlation must be due to chance, under the assumption that UHECR acceleration 
is not episodic in nature. 

9. SUMMARY AND CONCLUSIONS 

We have examined the 21 galaxies from the VCV catalog with z < 0.018 which are w i thin 3 .2° of the 27 highest- 
energy cosmic rays, and 8 additional AGNs from NED introduced bv IMoskalenko et alj (|2008f ). We find that 14 of 
the 21 VCV galaxies are AGNs and the other seven are either star forming or quiescent galaxies. X-ray or radio 
observations are needed to find out if any of these seven have obscured nuclear activity at an interesting level. Two of 
the 8 additional NED galaxies can be optically established as AGNs (one is a Seyfert 2, one a LINER, and 6 do not 
have adequate optical spectra to make a determination). The new optically identified AGN, WKK 4374, increases by 
one the number of UHECR- AGN correlations, because it is within 2.8° of a previously uncorrelated UHECR. 

We have determined the spectral type and bolometric luminosity of the 14 VCV AGNs and the newly identified 
AGN. Their luminosi ties range f rom 5 x 10 42 erg s _1 to 1 x 10 46 erg s _1 . Five are broad line, nine are narrow-line 
and one is a BL Lac. lHao et all (|2005h finds that the numbers of broad line and narrow-line AGNs are comparable at 
low and moderate luminosity, while at high luminosity, broad line AGNs far outnumber narrow-line ones. Thus the 
AGNs possibly producing the observed UHECRs are representative of the population of moderate luminosity AGNs 
both in luminosity and type. 

We have also identified and determined the bolometric luminosities of the Virgo AGNs and LINERs. Most have 
luminosities more than an order of magnitude below the lowest luminosity of any AG N correla t ed wit h a UHECR. (A 
large number of low luminosity AGNs and LINERs are known in Virgo, due to the iHo et all (|1995f) high-sensitivity 
spectroscopic survey of the near-by galactic nuclei in the northern hemisphere.) A minimum bolometric luminosity 
requirement for UHECR-accelerating systems may therefore explain the absence of UHECRs correlated with Virgo 
AGNs. If there is such a threshold, the significance of correlations between UHECRs observed by northern hemisphere 
observatories and VCV galaxies would be reduced in comparison to that observed by Auger South, due to a greater 
dilution of the VCV catalog in the northern hemisphere by large numbers of known very low luminosity AGN s. A 
threshold in the bolometric luminosity is consistent with the Giant AGN Flare scenario (|Farrar fc Gruzinovll2008f) . but 
is not predicted by it, since the properties of the remnant have not yet been modeled. 

Altogether, 21 of the 27 Auger UHECRs are within 3.2° of a Veron-Cetty Veron galaxy or the newly identified 
AGN. Of these, 17 of the 21 correlated UHECRs can be attributed to an identified AGN within 3.2° and z < 0.018. 
However few of the correlated AGNs satis fy the condition Lboi > 10 45 -E|o ergs -1 , required to confine a cosmic ray as 
it is accelerated (jFarrar fc Gruzin"ovi r2008) . If correlated AGNs with inadequate bolometric luminosities are in fact the 
sources of most of these UHECRs, their l uminosities presen t a pu zzle to the conventional picture of AGN acceleration 
in a continuous jet. The recent analysis of Moskal enko et al.l (|2008f) finds that few of the correlated AGNs have powerful 
radio jets or lobes, corroborating this conclusion. If the AGN-UHECR correlation is a real one, t hese results provide 
eviden ce in favor of the new mechanism of UHECR acceleration by giant AGN flares proposed by iFarrar fc Gruzinovl 
( 2008) , in which a modest AGN has an intense flare producing a jet of the required luminosity (initiated for instance 
by a stellar tidal disruption rapidly heating the accretion disk) and then reverts to a mild-mannered existence. 
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APPENDIX 
A. NOTES ON INDIVIDUAL OBJECTS 

The 21 VCV galaxies correlated with Auger UHECRs with energy above 57 EeV are comprised of one BL Lac 
(Centaurus A), 5 broad line AGNs (Seyfert 1 galaxies), 8 narrow-line AGNs (Seyfert 2 galaxies) and 7 (non-optical- 
AGN) emission line galaxies. In this section we give detailed information on the identification and properties of each 
of these galaxies. Work is underway to examine all of the VCV galaxies with z < 0.024, to determine their correct 
classifications and find the bolometric luminosities of those with active nuclei (IZ, GRF, N. Bursky-Tammam, E. Storm 
and A. Hoeft, in progress.) 

Centaurus A 

Centaurus A (NGC 5128) (|Israe1ll998ft is a well known, nearby (3.4 Mpc) BL Lac. Its nuclear acti vity is manifested 
by lar ge radio/X-ray jets - radio plumes extending to 250 kpc, and a compact circumnuclear disk. iRothschild et al.l 
(2006) report a 2-10 keV X-ray flux which varies between 1.69 - 3.23 x 10 -10 erg s _1 cm~ 2 in six observations between 
August, 1996 and February, 2004 with the PCU2 instrument on RXTE. 

Broad- Line AGNs (Seyfert 1 Galaxies) 

SDSS J03302-0532 (J033013. 26-053235. 9): The optical spectrum shows broadened Ha and H/3 lines. We fit 
m ulti-Gaussian mo dels to the spectrum. Galaxy subtraction is perfor med using the princip le compo nent analysis cod e 
of lHao et alj (120051) . and the narrow lines arc modeled as outlined in I Greene fc Hoi (|2004h ; see also lHo et ali (|1997d) . 
Following I Greene fc Hoi (|2005f ). the AGN l uminosity is measure d from the Ha line and the velocity dispersion is mea- 
sured from the Ha line width (for details see lGreene fe Hoi (|2007f l). The black hole masses for the SDSS objects are esti- 
mated from the broad Ha luminosity and line width using the relation M B h = 2.0x 10 6 ( 10 4g° s -i )°' 55 ( i(^™ H ° ) 2 ' O6 M0 
following [Greene fc Hoi (l2005h . We find: FWHM(Ha) = 5340 km/s, log(L(Ha)) = 40.83 erg s"\ log(L[OIII]) = 39.50 
erg s"\ M BH = 7.2 x 10 7 M , L bol /L Edd = 0.006. 

SDSS J03349-0548 (J033458. 00-054853.2): The optical spectrum shows broadened Ha and H/3 lines. We fit 
the spectrum as described for SDSS J03302-0532 and find: FWHM(Ha) = 4390 km/s, log(L(Ha)) = 40.10 erg s~\ 
log(L[OIII]) = 39.30 erg s" 1 , M BH = 6.6 x 10 7 M Q , L bol /L Edd = 0.015. 

Zw 374.029 (CGCG 374-029 ): A broadened Ha line can be seen in the optical spectrum with a broad Ha FWHM 
of 2048 km/s (jPietsch et al.ll2000f ). We have fit a single Gaussian to the narrow [OIII] line and find F[OIII] = 2.74 x 
10 -14 erg s _1 cm~ 2 . 

NGC 424: The broad Ha and H/3 lines were detected in polarized emission bv iMoran et~aT1 (|2000h . They report 
an Ha FWZI of approximately 12,000 km/s. However, since the broad lines are not detected in unpolarized emission, 
we use the narrow-line [OIII] flux, F[OIII] = 4267.0 x 10~ 16 erg s" 1 cm -2 (|Gu et al.ll2006h to derive the bolometric 
luminosity. 

4U 1344-60: iPiconcelli et alJ (|2006h report a broad Ha line with a FWHM of 4400 km/s. They classify it as an 
intermediate Seyfert (SI. 5) due to the absence of a broadened H/3 line. They report a de-reddened broad Ha flux of 
F(Ha) = 3.3 ± 0.8 x 10~ 12 erg s" 1 cm" 2 , and a 2-10 keV X-ray flux of iViokcV = 3.6 x 10" 11 erg s^ 1 cm~ 2 . 

Narrow- Line AGNs (Seyfert 2 Galaxies) 

NGC 5506: The optical spectrum for NGC 5506 has strong [0111] an d [Nil] narrow lines. The log( [OIII] /H/3) and 
log([NII]/Ha) line ratios , 0.88 and -0.09 respectively (|Kewlev et al.l l200D. place it firmly in the Seyfert 2 region of the 
BPT diagram. iGu et all (|2006h report an [OIII] flux, i^[OIII] = 1614.0 x 10" 16 erg s _1 cm" 2 and an X-ray luminosity, 
log(L2-iokcv) = 42.89 erg s" 1 . NGC 5506 also hosts circumnuclear H2O masers. 

IC 5135 (NGC 7130): The optical spectrum for IC 5135 has stro ng [OIII] and [NH ] narrow lines. The 
log( [OIII] /H/3) and log([NII]/ Ha) line rat i os, 6990 and -0.0231 respectively (jVaceli et al.lll997t ). place it in the Seyfert 
2 region of the BPT diagram. IShu et alJ (|2007h report an [OIII] luminosity, log(L[OIII]) = 42.55 erg s -1 . 

MRK 60 7 (NGC 1320): The log( [OIII] /H/3) and log([NII]/Ha) line ratios f or MRK 6 7, 1.0 458 and -0.1498 
respectively ([Vaceli et al.lll997t ) , place it in the Seyfert 2 region of the BPT diagram. IShu et all (|200l report an [OIII] 
luminosity, log(i[OIII]) = 41.05 erg s~ 4 . MRK 607 also hosts circumnuclear H2O masers. 

IC 4518A: We fit for the narrow lines in the optical spectrum of IC 4518A (courtesy of N. Masetti). The 
log( [OIII] /H/3) and log([NII]/Ha) line ratios (courtesy of J. Moustakas), 0.839 and -0.199 respectively, place it in 
the Seyfert 2 region of the BPT diagram. We fit a single Gaussian to the [OIII] line and obtain a flux of F[OIII] = 
1.20 x 10" 13 erg s _1 crrT~ 2 . 

NGC 1358: The log( [OIII] /H/3) and log([NH]/ Ha) line ratios 1.0542 and 0.3032 respectively (|Ho et al.lll997ah . 
place it in the Seyfert 2 region of the BPT diagram. IShu et all (|2007f ) re port an [OIII] lum inosity, log(L[OIII]) = 40.86 
erg s _1 . The SMBH mass is determined to be 6.56 ± 2.50 x 10 7 M (|Wu fc Hanll2001h . This implies an Eddington 
ratio of 0.029. 

ESQ 383-G18: T he log( [OIII] /H/3) and log([NII]/Ha) line ratios for ESO 3 83-G18, 9289 a nd -0.6327 respectively 
(|de Grijp et al.lll992f l. place it in the Seyfert 2 region of the BPT diagram. IGu et alJ (|2006l ) report an [OIII] flux, 
F[OIII] = 360.9 x 10" 16 erg s" 1 cm" 2 . 

NGC 4945: Optical emission from NGC 4945 is highly obscured. Even X-rays below 10 keV are absorbed by a 
column of Nh = 4.5 x 10 24 cm -2 , but NGC 4945 is one of the brightest Seyfert galaxies at 100 keV, indicating an 
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active nucleus (Madeiski et al. 2000). From 1-500 keV X-ray luminosity, Madeiski et al. (2000) estimate the bolometric 
luminosity to be ~2 x 10 43 erg s _1 . The mass of the SMBH has been determined from circumnuclear H 2O masers to 
be ~1.4 x 10 6 M Q , which implies that the Eddington ratio for NGC 4945 i s ~0.1 (|Greenhill et al.lll997t) . 

ESO 139-G12: The optical spectrum for ESO 139-G12 iMarquez et all (|2004D does not extend to [OIII] and H/3 
wavelengths, but the hint of a broadened Ha line indic ates that it is an AG N. Since the spectrum is not flux calibrated, 
we use the upper limit for the X-ray luminosity from iPolletta et all (|1996[ L L2-iokeV < 47.56 x 10 42 erg s _1 . 

Non- AGN Emission Line Galaxies 

NGC 2989: VCV identifies NGC 2989 as an HII galaxy. iGoncalves et all (|1999h (NB: Co-authors are Veron-Cetty 
and Veron) classify this as a pure HII galaxy, based on flux ratios, log( [OIII] /H/3) = -0.2304, log ( [Nil] / Ha) = -0.2757. 
This falls within but almost on the theoretical star formation line developed by iKewlev et al.l (j2001) . It may be a 
composite AGN-Starburst (transitioning) object with AGN contribution to the total luminosity < 20%. 

NGC 7591: The logj ' [QUI] /H/3) and log([NII]/Ha) line ratios, -0.345 5 and -0.3224 respe ctively 
(Moustakas & Kcnnicutt 2006), place it below the star formation line determined by IKauffmann ~et al.l (120031 ) . The 
infrared luminosity, Log(L/_R) = 11.12 L , indicates tha t it is a luminous infrared galaxy (LIRG). ISturm et aH(|2006D 



classifies it to be a Starburst galaxy while IVeilleux et"aTI (|1995l ) classifies it as a LINER, resulting from shocks and not 
from AGN accretion activity. 

NGC 1 204: The (de-redden ed) log( [OIII] /H/3) an d log([NII1/Ha) line rat ios fo r NGC 1204 - 0.47 and -0.11 re- 
spectively (IVeilleux et al.lll995l). pla ce it between the Kauffm ann et al.l (|2003f ) and IKewlev et al.l (|2001f ) lines in the 
LINER re gion. IVeilleux et al.l (fl995l ) give reasons to believe that the emission l ines are a resu l t of s hocks instead of 
accretion. ISturm et al.1 (|2006l ) also classify it as a star forming galaxy. However. ICorbett et al.1 (|2003f ) . based on their 
radio and spectroscopic COLA survey, claim this galaxy lies on the border between Seyfert and LINER and classify it 
as a Seyfert even though no radio core was detected and the radio continuum is dominated by star formation. They 
place a limit on the AGN component, of 25% of the total e mission of the gala xy. 

NGC 5244: While the optic al spectrum f o r NG C 5244 (jMoran et al.lll996h has [Nil] and Ha emission lines, [OIII] 
and H/3 lines are not apparent. [Moran et alj (|1996f ) classifies it as an HII galaxy and VCV also identifies it as an HII 
galaxy. 

IC 5169: (|Kewlev et al.ll2001ft classifies this as an HII galaxy based on the following line ratios, log( [OIII] /H/3) 
= 0.20, log([NII]/Ha) = -0.1 70, log([OIl63 0/Ha ) = -1.49, log([SII]6717,31/Ha) = -0.55. While it falls between the 
IKauffmann et al.1 (|2003f) and IKewlev et all (|200lD lines in the log( [OIII] /H/3) vs. log([NII]/Ha) diagram, it falls well 
within the star formation line in log( [OIII] /H/3) vs. log([SII]/Ha) and log( [OIII] /H/3) vs. log([OI]/Ha) diagrams, 
confirming that it is not an o p tical A GN. 

NGC 7135: Uoguet et all (|2001h classify NGC 7135 as a star forming galaxy based on the fact that the [OIII] 
emission line is very weak, EWQOIII]) = -0.7A, and that strong stellar absorption features are visible, determining it 
to be "an early-type galaxy with no active nucleus" . 

Q 2207+0122 (PC 2207+0122): Although the optical spectrum for Q 2207+0122 (courtesy of D. Schneider) 
shows emission lines, Q 2207+0122 is just an emission line galaxy and not an AGN (confirmed by D. Schneider, 
private communications). The log( [OIII] /H/3) and l og([NII]/Ha) line ratios (courtesy of J. Moustakas), 0.815 and 
-1.581 respectively, place in the star formation region IKauffma nn et al.l (|2003l ). This galaxy is a low metallicity dwarf 
galaxy (J. Moustakas, private communication). 

B. DETAILS ON THE AGNS IN THE VIRGO CLUSTER 

The details for the VCV galaxies in the Virgo Cluster are lis ted in the table belo w, in order of distance from the 
center of the cluster at a J20 oo = 12:26:32.1, 6j 2000 = +12:43:24 (|Binggeli et al.lll987l ). The [OIII] and 2-10 keV X-ray 
luminosities, as well as the bolometric luminosities derived from them, are given for the AGNs not listed as HII galaxies 
in VCV. Note that the bolometric luminosities derived from [OIII] luminosities are systematically higher than those 
derived from 2-10 keV X-ray luminosities. There is some evid e nce th at the [OIII] to bolometric luminosity conversion 
depends on luminosity or Eddington ratio (e.g., iNetzer et al.l (l2006[)). S ince optical emission line fluxes, particularly 
Ha, may be contaminated by non-nuclear sources re.g. JHo et alJ (|2001h ) we have chosen to use the hard X-ray flux 
as a bolometric indicator whenever possible. 

C. NOTES ON THE AGNS FROM THE NASA EXTRAGALACTIC DATABASE 

iMoskalenko et all (|2008l ) lists eight additional AGNs within 3.2° of the Auger CRs with z < 0.018, based on the 
classifications in the NASA Extragalactic Database. Some of them were determined to be AGNs after the publications 
of the VCV 12th Edition. Others were selected as AGNs mainly via X-ray, infrared, or radio criteria and show only 
some evidence of nuclear activity in the optical, which may explain their absence from VCV. One of these AGNs, WKK 
4374, correlates with a CR (Auger year = 2007, Auger day = 186; E=64 EeV), which was previously uncorrelated. 
Furthermore, two of the additional AGNs, NGC 7626 and NGC 2907, correlate with two CRs, (2006, 185; E=83 
EeV) and (2007, 84; E=64 EeV), respectively, whose VCV counterparts were shown in this paper to be non-nuclear 
dominated emission line galaxies. We perform a literature search on the extra AGNs and list their optical properties 
below. 

WKK 4374 (IGR J14515-5542): From optical spectra of INTEGRAL sources. iMasetti et all (|2006l ) determined 
WKK 4374 to be a Seyfert 2 galaxy. They report an absorption-corrected [OIII] flux of 6.4 ± 1.6 x 10~ 14 erg s -1 cm~ 2 . 
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TABLE 4 
Details on the Virgo AGNs in VCV 



Name 


K A 


Dec 


z 


VC V Olass 


Hol997Class 


Log(L[OHIJ) 


Log(L 2 -10kev) 


LbolAOlII] 


Lbol, 2 -lOkeV 




(J2000) 


(J2000) 








erg s — 1 


erg s — 1 


erj 




— 1 


erg s — 1 


NGC 4388 


186.44 


+12.66 


0.008 


Slh 


SI. 9 


40.35 


41.72 a 


7.85 


X 


10 43 


1.05 x 10 48 


NGC 4374 


186.27 


+12.89 


0.003 


S2 


L2 


38.45 


39. 12 6 


9.78 


X 


10 41 


2.64 x 10 40 


NGC 4438 


186.94 


+13.01 


0.004 


S3b 


LI. 9 


38.80 


39.22 c 


2.22 


X 


10 42 


3.32 x 10 40 


NGC 4486 


187.71 


+12.39 


0.004 


S3 


L2 


39.07 


40.75 c 


4.15 


X 


10 42 


1.12 x 10 42 


NGC 4477 


187.51 


+13.64 


0.005 


S2 


S2 


38.82 


39.65 a 


2.30 


X 


10 42 


8.93 x 10 40 


NGC 4501 


188.00 


+14.42 


0.007 


S2 


S2 


39.10 


39.59" 


4.38 


X 


10 42 


7.78 x 10 40 


NGC 4552 


188.92 


+12.56 


0.001 


S2 


T2 


38.05 


39.55 d 


3.94 


X 


10 41 


7.10 x 10 40 


NGC 4550 


188.88 


+12.22 


0.001 


S3 


L2 


38.15 


< 37.68 d 


4.95 


X 


10 41 


< 9.57 x 10 38 


NGC 4569 


189.21 


+13.16 


0.004 


S 


T2 


39.27 


39. \\ h 


6.54 


X 


10 42 


5.2 x 10 4U 


NGC 4548 


188.86 


+14.51 


0.002 


S3 


L2 


38.11-'' 




4.54 


X 


10 41 




NGC 4579 


189.43 


+11.82 


0.005 


S3b 


S1.9/L1.9 


39.42 


41.15 c 


9.25 


X 


10 42 


2.83 x 10 42 


NGC 4168 


183.07 


+13.21 


0.008 


SI. 9 


SI. 9 


37.91 


39.87 a 


2.84 


X 


10 41 


1.48 x 10 41 


NGC 4383 


186.36 


+16.47 


0.005 


H2 
















NGC 4192 


183.45 


+14.90 


0.004 


S3 


T2 


38.12 


39.58 e 


4.57 


X 


10 41 


7.6 x 10 40 


NGC 4639 


190.72 


+13.26 


0.001 


S1.0 


S1.0 


38.41 


40.22 a 


8.96 


X 


10 41 


3.32 x 10 41 


NGC 4450 


187.12 


+17.08 


0.006 


S3b 


LI. 9 


38.78 


40.35 e 


2.10 


X 


10 42 


4.46 x 10 41 


NGC 4472 


187.45 


+8.00 


0.003 


S2 


S2 


37.81 


< 39.32 a 


2.25 


X 


10 41 


< 4.18 x 10 40 


NGC 4694 


192.06 


+10.98 


0.004 


H2 


HII 














IC 3576 


189.16 


+6.62 


0.003 


H2 
















NGC 4698 


192.10 


+8.49 


0.003 


S2 


S2 


38.81' 




2.25 


X 


10 42 




NGC 4303 


185.48 


+4.47 


0.005 


S2 


HII/L 9 


39.13 


39. 16 9 


4.72 


X 


10 42 


2.89 x 10 40 


NGC 4412 


186.65 


+3.97 


0.008 


S2 


S2 h 














NGC 3976 


178.99 


+6.75 


0.008 


S2 


S2 


39.12 




4.61 


X 


10 42 




NGC 4636 


190.71 


+2.69 


0.003 


S3b 


LI. 9 


38.10 


39.25 c 


4.41 


X 


10 41 


3.56 x 10 40 


MRK 52 


186.43 


+0.57 


0.007 


H2 
















NGC 4772 


193.37 


+2.17 


0.003 


S3b 


LI. 9 


38.45 




9.86 


X 


10 41 




SDSS J12122+0004 


183.06 


+0.07 


0.008 


SI 


HII 














(MRK 1313) 
























NGC 4418 


186.73 


-0.88 


0.007 


S2 


S2 h 




39.25 i 








3.56 x 10 40 


SDSS J12580+0134 


194.51 


+1.58 


0.004 


S2 


HII 


38.19 




5.45 


X 


10 41 




(NGC 4845) 
























MRK 1308 


178.55 


+0.14 


0.004 


H2 

















Note. — S=Seyfert, L=Li ner, T=Transitio n object. The classifications in column 6 and the [OIII] luminosities in column 7 (except 
where marked) are taken from Ho ct al. (1997a). 

References. — (a ) IPanessa et~all [[20061) (b) IHo et al.1 [120011), ( c) IGonzalez-Martm et all J200^>. (d) ISatvapal et al.l (t2005Tl . (e) 
ITerashima et al.l [12003) . (f) IShields et al] 1120071) . (g) IHorst et al.1 H20081) . (h) IKewlev et al.l H200lT) . (i) IKandalvanl [1200511 

This gives an [OIII] luminosity of 5.4 x 10 40 erg s _1 and implies a bolometric luminosity of 1.9 x 10 44 , within the range 
found for the correlated AGNs from VCV; A bo i = 0.46. 

NGC 5140 and N GC 2907: These two galaxies were determined to have narrow LINER-like emission by 
iMauch &: Sadlerl (|2007f ) from visual inspections of the 6dFGS optical spectra. No emission line fluxes or ratios have 
been published. 

NGC 7626: This is a radio loud galaxy. IHo et all (|1997a[ ) determine it to have a possible L2 nucleus (the classifi- 
cation is characterized as highly uncertain), with an [OIII] luminosity of 3.3 x 10 38 erg s _1 (100% uncertainty), which 
implies a bolometric luminosity of 1.2 x 10 42 erg s — 1 , giving Aboi = 0.0017. Since some correlated galaxies are chance 
coincidences, this may be one of them. 

PKS 1343-60 (Centaurus B): Cen B is a well-known FRI radio galaxy. However, while an optical counter-part 
for the active nucleus has been detected (Wes t fc Tarenghilll989f ). the optical emission is highly absorbed and no clear 
emission lines o r line ratios cou l d be o btained from the optical spectrum beyond the detection of the Ha- [Nil] complex. 

NGC 5064: iBonatto et all (jl989l ) report a [Nil] /Ha value of 1.051 and [SII]/Ha value of 0.461. No [Oil], H/3, or 
[OI] line fluxes were reported. This indicates that NGC 5064 is a LINER. The authors determine this LINER to be 
nuclear powered rather than shock powered. 

IRAS 13028-4909: This galaxy is from a set of X-ray and IR selected AGNs (jKirhakos fc Steinerl fl990t ) . An 
[Nil] /Ha ratio of 1.28 was reported. No other line fluxes or ratios wer e given. 

AM 1754-634 NED03: Shows a bright optical nucleus. However [Fairalll (|1981l) reports that the 5007 and 4954 
[OIII] lines were of equal strength (instead of in a 3:1 ratio), Ha was barely detected, and no other emission lines were 
detected. This casts doubt on the strength of its nuclear activity based on optical data. Data from other wavelengths 
has not been published. 
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